Emergent Universe Scenario in the 
Einstein-Gauss-Bonnet Gravity with Dilaton 

B. C. Paul *and S. Ghose t 
Physics Department, North Bengal University 
Dist : Darjeehng, PIN : 734 430, India. 



Abstract 

We obtain cosmological solutions which admit emergent universe 
(EU) scenario in the framework of Einstein Gauss-Bonnet gravity cou- 
pled with a dilaton field in four dimensions. The coupling parameter of 
the Gauss-Bonnet terms and the dilaton in the theory are determined 
for obtaining an EU scenario. The corresponding dilaton potential 
which admits such scenario is determined. It is found that the Gauss- 
Bonnet (GB) terms coupled with a dilaton field plays an important 
role in describing the dynamics of the evolution of the early as well 
as the late universe. We note an interesting case where the GB term 
dominates initially in the asymptotic past regime, subsequently it de- 
creases and thereafter its contribution in determining the dynamics of 
the evolution dominates once again. We note that the Einsteins static 
universe solution permitted here is unstable which the asymptotic EU 
might follow. We also compare our EU model with supernova data. 
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1 Introduction: 



Recent cosmological observations predict that the universe is passing 
through an accelerated phase of expansion [1]. Although the Ein- 
steins General theory of Relativity (GTR) with usual matter fields 
of the standard model of particle physics admits early inflation fairly 
well [21 [3], it fails to accommodate the recent accelerating phase. This 
led in modern cosmology to explore a suitable theory which accom- 
modates observational prediction. Perhaps the present accelerating 
phase of expansion may be realized either modifying the matter sec- 
tor or the gravitational sector including higher derivative terms that 
are relevant in the low curvature limit. It is generally accepted that 
inflation is an essential ingredient to construct cosmological models in 
modern cosmology. The early inflation can be realized consistently in a 
semi-classical theory of gravity |3l i4j . It may be mentioned here that 
Starobinsky [5] obtained inflationary solution in a higher derivative 
theory of gravity long before the advent of inflation actually realized. 
The gravitational action of the Starobinsky model corresponds to a 
theory which contains curvature squared terms in the the Einstein- 
Hilbert action. However, the efficacy of inflation is known only after 
the seminal work of Guth [3] who used phase transition mechanism to 
obtain such scenario in order to resolve some of the outstanding prob- 
lems of Big bang model. The idea of inflation is very much attractive 
which has been implemented in various theories leading to a number 
of versions of inflation in the literature. However, the late accelerating 
phase is very recently predicted from observational data and a suit- 
able explanation is yet to come out. To address the issue a number 
of proposals came up (i) with a suitable modiflcation of the matter 
sector by incorporating exotic kind of matter flelds [6] or (ii) with a 
modification in the gravitational sector by adding curvature squared 
or its inverse power terms that are effective at low curvature limit [7] 
to the Einstein-Hilbert action. It has been shown [7] that the present 

4 

accelerating phase of the universe may be obtained by adding a ^ 
term to the Einstein-Hilbert action. However, Einstein-Hilbert action 
with the inverse curvature correction term only is not enough as it is 
not free from shortcomings. Later, it has been shown that the above 
shortcomings of the theory may be removed by adding a polynomial in 
Ricci scalar i.e., curvature squared term to the proposed gravitational 
action [8j. It was shown long ago by Zweibach [9] that the string cor- 
rections due to Einstein action up to first order in the slope parameter 
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and fourth power of momenta should be proportional to Gauss-Bonnet 
(GB) terms (where GB = R^^^sR^''^^ - ^R^,uR^'' + R^)- However, 
it was realized subsequently that the field redefinition theorem of 't 
Hooft and Veltman [10] may be applicable in this case. On the Ein- 
stein's shell {R^u = 0), an action with curvature squared term of the 
form R + aR?^^, + hR^ may be transformed into R itself (neglecting 
higher order terms) by field redefinition : 

, P a + 2h 

9fj.u — 9fiu + aRfj^u + 9fiu ^_ j^ R 

where D represents the number of dimensions. Subsequently Deser 
and coworkers pTj have shown that on the linearized Einstein shell, 
the actions R+a' {GB) and R+a' R^^^^^^ (here a' is the inverse of string 
tension) are not different and this result generalizes to all higher-order 
ghost terms. GB terms arise naturally as the leading order of the a' 
expansion of heterotic superstring theory [12]. It is known that GB- 
terms in the higher dimensions leads to ghost free propagator. It may 
be mentioned here that in the framework of 4 dimensions, the GB 
terms do not contribute in the dynamics of evolution. However, GB 
terms if coupled with a dilaton field in the action then the combina- 
tion plays an important role in the dynamics of the evolution via the 
dilaton field. It has been shown also that GB combinations coupled 
with scalar field plays an important role for avoidance of singularity 
in a string induced gravity jl3j . Later, in the braneworld scenario, 
it has been shown [H] that the naked singularities may not occur if 
dilaton with a Gauss-Bonnet term are considered. The issues of fine 
tunning in a theory with a dilaton field and GB-terms interaction are 
also taken into account to study in details in Ref. [141 [T5] . Recently, 
cosmological models with dark energy of the universe are probed in 
the Einstein-Hilbert action with GB term and it is known that the 
theory accommodates the new form of energy [16]. It has been shown 
|17j that scalar field that enters into the coupling of GB terms in the 
action plays an important role which may be used to explain different 
phases of expansion of the universe including the present accelerating 
phase. It is also shown that GB terms with dilaton admits accelerating 
cosmologies in the framework of higher dimensions |19j . The Gauss- 
Bonnet terms in the Einstein-Hilbert action are used to obtain new 
black holes solution [18] and Kaluza-Klein space-times |2H [22]. Paul 
and Mukherjee |24j earlier noted that a Gauss-Bonnet term in higher 
dimensions leads to 3 a 4-dimensional universe at a later epoch with 
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many good features considering the sign of the couphng parameter 
different from that one usually gets from the low energy limit of string 
theory. The model also gives a satisfactory explanation of the small- 
ness of the effective four dimensional cosmological constant. Recently 
Einstein-Hilbert action with a combination of higher order curvature 
terms e.g., GB terms including dilaton are employed to study the 
present acceleration of the universe |171ll8j. Therefore, the GB-theory 
has a rich structure that needs to be explored. In this paper we ex- 
plore emergent universe scenario in the Einstein-Hilbert action with 
GB terms coupled with dilaton field in the frame work of four dimen- 
sions. We employ the GB term with dilaton in four dimensions to 
obtain an emergent universe scenario [241 [25] . 

Earlier Harrison [26] obtained a cosmological solution with radia- 
tion in the presence of a cosmological constant in a closed model of 
the universe which asymptotically approaches to Einstein static uni- 
verse but the scenario does not exit inflationary phase. Recently, Ellis 
and Maartens [27] obtained similar cosmological solution considering 
a minimally coupled scalar field with a special choice of its poten- 
tial where the universe exits from its inflationary phase followed by 
reheating. Subsequently it was shown by Ellis et al. [23j that the po- 
tential required to obtain such scenario may be obtained naturally by 
a conformal transformation of Einstein-Hilbert action with i?^-term 
for a proper choice of its coupling constant. The model incorporates 
an asymptotically Einstein static universe in the past and it evolves 
to an accelerating universe in the framework of a closed model of the 
universe. This model is usually known as emergent universe model. 
The salient features of an emergent universe scenario is that there is 
no time like singularity, it is ever existing and it approaches a static 
universe in the infinite past (t — > — oo). It is interesting to construct 
an emergent universe model as it is capable of solving some concep- 
tual issues of the standard Big bang model. The asymptotic Einstein 
static universe at some stage enters into the standard Big bang phase 
and might have features precisely known to us. The possibilities of 
an emergent universe scenario have been studied recently in a number 
of theories [HIIISIIITIISHIIMIIMIEIIEIIMIIM] because it permits 
a universe which is ever existing and large enough so that the space- 
time may be treated as classical entities. Recently, Mukherjee et al. 
[24] obtained an emergent universe (EU) scenario in a flat universe in 
the modified Starobinsky model. Subsequently, Mukherjee et al. [25] 
proposed a general framework for such an EU scenario in GR with 
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a mixture of matter and exotic kinds of matter that are prescribed 
by an equation of state (EOS) : p = Ap — B^fp^ where A and B are 
constants. The EU scenario can be reahzed for some possible soup 
comprising primordial compositions of matter that are permitted by 
the EOS. It admits existence of exotic matter and dark energy in ad- 
dition to radiation/dust Since recent cosmological observations 
indicate that our universe is almost flat, the above emergent model 
of the universe is explored in a flat universe context. Subsequently 
the EU scenario proposed by Mukherjee et al. [Ml EH] are explored 
in the context of various theories and found that the scenario can be 
realized fairly well |3 H 132^ 133]. The purpose of the paper is to examine 
EU scenario in a gravitational action with GB-terms coupled with a 
dilaton field. The plan of the paper is as follows: in sec. 2, we present 
the gravitational action and set up the relevant field equations , in sec. 
3, cosmological solutions with emergent universe scenario are derived. 
Finally, in sec. 4, we summarize the results obtained. 



2 Action and the Field Equations: 

We consider a gravitational action given by 

/ = - / + fmR^.u^sR^''^' - ^RapR""" + R^) + d^x 

(1) 

where the Greek indices /U, u represents (0, 1, 2, 3), /(0) represents 
the coupling factor of the Gauss-Bonnet (GB) term and GB is a com- 
bination of squared terms of Riemann tensor, Ricci tensor and Ricci 
scalar ( GB = R^^^sR^'"^^ - ^R^pR"''^ + R^), 9 represents the 4 di- 
mensional metric, SttG = k? and represents the Lagrangian for the 
dilaton field. The corresponding Lagrangian for the dilaton field is 
given by 

= -m d^^d^c^ - (2) 

where ^(i;^) represents the coupling parameter for the field in the grav- 
itational action and l^((/>) represents potential of the dilaton field. 

We consider a flat, homogeneous and isotropic Robertson- Walker 
(RW) metric with scale factor a(t), which is given by 
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The action given in (1) with the RW metric (3) yields the following 
field equations : 



C(0)</)2 + y('/')-24/'(<^)A 



(4) 



2-+ - 



(5) 

once again varying the action (2.1) with respect to the dilaton field 0, 
we get 



a 2 ^ 2^ 



a^a 



12/'(</>)^, 



(6) 



where the over dot implies derivative w.r.t. time and prime (') repre- 
sents differentiation with respect to the field (p. In the above, out of 
the three eqs. (4)-(6), only two are independent, as eq. (2.6) can be 
derived from eqs. (4) and (5). It is evident that there are altogether 
five unknowns namely, a{t), (f), V{(j)), ^{4>) and /(</>) in the above field 
equations. In order to solve the equations additional assumptions are 
can be made. Let us first assume that 



V, 



(7) 



where is a constant. The above assumption leads to a relation 

f{(j)) = rjt{(f)) + 77o, where r/„ is a constant. The coupling parameter 
/(0), therefore, grows with time. Consequently the effect of the GB 
terms becomes more and more important at late time, which may be 
useful for describing the present acceleration of the universe. Using 
the constraint given by eq. (7) in eqs. (4) and (5), one gets 



2H + 3iJ^ = -K^ 



4(0)(/)2 + y(0) - 247jH^ 



(8) 
(9) 



where H = ^ is the Hubble parameter. Now eliminating V{(f)) from 
eqs. (8) and (9), we get 



^(0)(/>2 + SrjHH - 4:7] 



0. 



(10) 



6 



The dilaton potential can be obtained from eqs. (8) and (9) eliminat- 
ing ^((j)) (fP' , the corresponding potential becomes a function of Hubble 
parameter which is given by 

V{H) = ^H'^ + 2{)i]H^ + \{l + SriK^H)H, (11) 

the above potential may be expressed as a function of </> once the Hub- 
ble parameter {H) is known in terms of (j), which will be determined 
in the next section. We use eqs. (8)- (9) to determine the coupling pa- 
rameters ^{(j)) and /((/)) in the theory for an EU cosmological solution 
following Ref. [17]. The set of eqs. (lO)-(ll) contain four unknowns, 
therefore, two more ad hoc assumptions can be made to obtain a con- 
sistent cosmological solution. In the next section we begin with a 
known scale factor, which permits an emergent universe scenario. 



3 Cosmological Solutions 



Let us consider the evolution of the scale factor of the universe in the 
form 



a{t) 



A + e 



at 



(12) 



where Oq, a, (3 and A are positive constants. It gives an EU scenario 
as have been obtained in Refs. |241 125j. The Hubble parameter cor- 
responding to eq.(12) satisfies a first order differential equation given 
by 

H = aH - (3H^. (13) 

From eqs. (8)- (10), the dilaton coupling and the dilaton potential are 
now can be determined. The field equations are highly non-linear. 
Consequently it is not simple to obtain a general form of ^{4>) and 
V{(f)) in terms of the dilaton field. However, those parameters may 
be determined in terms of the Hubble parameter following Ref. |17j . 
which are 



C{H) = UrjK^ (1 + 213) H^+(p- 8K^ar]) - aH 

^2^2 L V / . 

V{H) = —H"^ + 2{)r]H^ + — (l + 8K%i/) iaH - pH"^ 



(14) 



(15) 



For simplicity, as a special case let us consider f3 = 8r]K?a in the 
above. the corresponding dilaton field potential can be expressed in 
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terms of the field. There is another freedom to assume here as the 
number of unknowns are one more than the number of relevant equa- 
tions. Therefore, we look for emergent universe scenario for different 
behaviors of the dilaton field in the next subsection. 

3.1 Case I : 

For an increasing dilaton field (j) = 4'o&°'^, with the GB coupling terms 
f{(j)) = ^ In (j) , dilaton coupling is given by 

We note the following: (i) ^(0) oo when d — > i.e. — > 
and (ii) ^(0) ^ at two points (a) for Hi = ^ and (b) for H2 = 

/3(i+2/3) • P^°^ variation of the dilaton coupling (^) with as 
shown in the fig.l. It is evident that initially the dilaton coupling 
begins with negative value (phantom like property) but in course of 
its evolution (^{(j)) becomes positive and almost constant after attaining 
a peak. Thus in this model to begin with one can start with a field 
having negative kinetic energy, an interesting field which behaves like 
phantom [35], now-a-days it is considered as one of the candidate of 
dark matter and avoids the singularity. The potential is given by 



The Hubble parameter in this case is related to field as 



We determine A = 1 , (f>o = and the corresponding evolution of the 
scale factor becomes 



Oo 1 + e 



1 

at /3 



(19) 



3.2 Case II 



For a decreasing dilaton field <j) = 4>oe with the GB coupling term 
/((/)) = —^ln<j), the dilaton coupling is given by 
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Figure 1: ^{(p) Vs. <p curve for k = 1 and (3 = 2 



In fig. 2 we plot variation of dilaton coupling ,^(0) with It is 
evident that ^{(p) begins with a negative value then attains a minimum 
thereafter it increases as (j) decreases. The potential for the dilaton 
field is obtained as 



a 
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y /? + {2(3 + 3)0 + (t? 



k2(/3 + 0)3 

The Hubble parameter is given in terms of the dilaton as 

a 



H 



the corresponding scale factor is 



(3 + ct> 



(21) 



(22) 



a = an 



4>o + Pe 



at 



(23) 



3.3 Case III : 



For a slowly varying field (j) = ^ In t, with /((/>) = r] e°"^, similar to that 
one expects in the string theory framework [34], the dilaton coupling 
is given by 

2 



e('A) = ^ 



2a 



aH 



In- 



H 



and the dilaton potential is 



V{H) 



1 



aH + + 



13 (h 



a 



a-PH 



13] 



(24) 



(25) 




Figure 2: ^(0) vs. (j) curve for k = 1 and (3 = 2 



where the Hubble parameter (H) is given by 



ae 

H = -. 26 

1 + pe'"''"^ 

Fig. 3 shows variation of .^(0) with <j) for (3 = 5. Here the dilaton 
couphng decreases as the dilaton increases to begin with from a nega- 
tive value initially, attains a minimum, thereafter, it increases sharply. 
GB terms become important at late time in this case. The potential 
becomes flat as (( ^ oo which is shown in fig. 4. 



3.4 Case IV : 

For a dilaton field (p = H, considering variation of the coupling of GB 
term as f = ijt. We obtain : 



with dilaton potential given by 



/3(1 + 2/?) 



2a 



a 



V{<l>) 



/3(5-2/3) . 



.a(f) + 3(/>^ + 
^ \ 2a 



(27) 



(28) 



The fig. 5 shows the variation of S,{(p) vrs. 0, for /3 = 5. We note 



that the coupling parameter becomes undetermined when 



/3- 



In 



this case it is necessary to begin with an initial field which greater 
than the above limiting value. In that case the coupling parameter 
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;ure 4: V{(j)) vs. (j) curve for a = 0.5 and (3 



always remains positive definite and it decreases from a large value 
to zero. Here the GB terms do not contribute at late time as the 
coupling ^ — > 0. The corresponding dilaton potential required for the 
EU model is shown in fig. 6. It has a minimum which is negative 
definite. The dilaton potential for /3 = |, becomes 



In this case the universe evolves as 



2 

5 i \ 5 



a = ao(l + -e"*) . (30) 



3.5 Case V : 

We now consider a special case /((/)) = r?o + 'n4'^ to obtain EU sce- 
nario. In this case the dilaton evolve as = ±-v/t with the coupling 
parameter 

m = K + + 2/3)^') • (31) 



The dilaton field potential is 

V{<P) = ^,{c.H + + H') (32) 

where i?=^^0^. 

Fig. 7 shows the variation of the dilaton coupling parameter with 
field, which is interesting. The variation of the potential is also shown 
in fig. 8, which has two flat regions one in the early era and the 
other in the late era respectively. This is new and interesting as it 
permits both early inflation and late acceleration. In Ref. [23| it 
was shown that a particle creation may occur during a phase when 
the Hubble parameter varies slowly. In this case the GB coupling 
parameter /(</>) — > ±oo as t — > ±00. The GB combination might have 
dominated in the early epoch of evolution which eventually decreases 
at later epoch (for rjo = 0) corresponding to a minimum of f{(j)) (say 
at t = 0) thereafter it increases which is shown in fig. 9. Thus GB 
terms play an important role in the early era which is subsequently 
important once again at late era contributing to the dark energy [8]. 
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Figure 5: ^{(/)) vs. (/) curve for a = 0.5 and (3 
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Figure 6: V{(j)) vs. (j) curve for a = 0.5 and (3 
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Figure 9: /(0) vs. curve for rjo = and 77 = 5 



14 



4 Stability of Einstein's Static Universe 



The EU scenario is characterized by the scale factor a{t) = (oq + 
P^at-ji/is g^ygj^ \yi (12). It may be pointed out here that the universe 
(i) begins with singularity if /? < 0, (ii) begins with singularity and 
asymptotically approaches Einstein static (ES) universe at late time if 
a < and P ^ 0, and (iii) spends infinite time near the Einstein static 
universe but pulls away and ends in an infinite inflating epoch if a > 
and P > 0. It is pointed out that the first two cases lead to unstable 
solution [37]. It may be pointed out here that: (i) The existence of 
the ES universe in fourth order theories of gravity and stability of the 
ES has been studied in |39j. (ii) Stability of ES in f{R) gravity has 
also been considered in [4^. (iii) It was also noted that ES universes 
are unstable in generic f{R) models [41]. To analyze the stability of 
Einsteins Static (ES) universe (which is an asymptotic past solution) 
[38j in the theory we begin with a pair of differential equations in a 
and H which are given by 

d = aH. (33) 

H = aH-l3H'^, (34) 

The above equations form an autonomous system which can be ana- 
lyzed by the standard technique. The Einstein static universe solution 
corresponds to the critical point of the system (oo,0). The ES uni- 
verse in this case is unstable for a > 0. It may be pointed out here 
that the stability of the Einstein Static universe under inhomogeneous 
perturbations has recently been studied in Refs. |42j . 



5 Distance modulus curve: 

We now probe late universe in the Emergent Universe scenario taking 
into account the observational results available from supernova. The 
distance modulus is ^ = Slogd^ -|- 25 where di is the luminosity 
distance (in the unit of mega parsecs), given by 

dL = ri{l + z)a{tQ) (35) 

where ri is given by : 

dr r^o rlf 

ar I at ^^g^ 



Vl - kr'^) Jti 
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Figure 10: /u vs. z curve for Emergent universe and supernova data. 



We consider the scale factor for a EU scenario as was given in (12), 

a(t) = ao(cT + e"*r, (37) 

where ao,o" and lo are constants. At late time the exponential term 
dominates and one can write 

a(t) ~ aoe"'^* (38) 

Since we are considering a flat universe, A; = 0, eq. (36) yields 

n = / ^TT (39) 
hi a[t) 

Using the scale factor (37) one obtains an expression for dL which is 
given by 

(40) 

Note that the final expression for di^ does not depend on ao,cr and 
ijj. It is now possible to determine [i{z) numerically for an Emergent 
Universe scenario at different values of redshift parameters (z). The 
observed values of //(z) at different z parameters p] along with that 
obtained from the present theory are given in Table : 1. 
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z 


Supernova /x 


EU u 


0.038 


36.67 


36.0438240 


0.014 


33.73 


33.7609468 


0.026 


35.62 


35.1945228 


0.036 


36.39 


35.9222306 


0.040 


36.38 


36.1593860 


0.050 


37.08 


36.6647158 


0.063 


37.67 


37.1932883 


0.079 


37.94 


37.7171620 


0.088 


38.07 


37.9694771 


0.101 


38.73 


38.2944627 


0.160 


39.08 


39.4068091 


0.240 


40.68 


40.4320839 


0.300 


41.01 


41.0192423 


0.380 


42.02 


41.6622327 


0.430 


42.33 


42.0079418 


0.490 


42.58 


42.3808310 


0.526 


42.56 


42.5866208 


0.581 


42.63 


42.8794593 


0.657 


43.27 


43.2483587 


0.740 


43.35 


43.6128241 


0.778 


43.81 


43.7684760 


0.828 


43.61 


43.9639519 


0.886 


42.91 


44.1787963 


0.949 


43.99 


44.3993095 


0.970 


44.13 


44.4701090 


1.056 


44.25 


44.7473544 


1.190 


44.19 


45.1438747 


1.755 


45.53 


46.4859129 



6 Discussion 

We obtain emergent universe scenario in a modified theory of gravity 
with Gauss-Bonnet term coupled to a dilaton field in four dimensions. 
We look for an emergent universe (EU) scenario here in a spatially 
flat universe. In the EU model, the universe was in an almost static 
state in an infinite past which eventually evolves into an inflationary 
stage later on. The Einstein static universe permitted here is found 
unstable. We note that emergent universe scenario obtained in Refs. 
\14:\ [25] can be implemented in a modified theory of gravity with GB 
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terms quite successfully with a suitable dilaton and GB coupling pa- 
rameters. The parameters are determined here and the corresponding 
models are presented in cases I to V in section 3. The result ob- 
tained here once again supports the view that the GB combination 
might have played a crucial role for driving both early inflation and 
late acceleration. A new solution is noted here in which the coupling 
parameter of the GB terms becomes important in the early and late 
era. This is an interesting solution with a new dilaton potential which 
is shown in fig. 8. The late acceleration of the universe may be ex- 
plained in this framework quite successfully as is shown in fig. 10 by 
the distance modulus curve. We also note cosmological solution where 
the dilaton field behaves like phantom |35j in the early era but it tran- 
sits to non-minimally coupled scalar field in the later epoch in all the 
cases except in one case, where (p = H and /(</>) increases linearly with 
time. However it is evident from fig. 6 that there is a regim where 
V{(f)) remains negative. We compare our Emergent Universe model 
recent SNela data [Ij at late universe. Fig. 10 show a comparative 
study of ^i{z) vs. z curve obtained from observation and that obtained 
theoretically from Emergent universe model. It may be pointed out 
here that between an exponential and accelerating phase, there might 
exist a phase of particle creation era which is pointed out in Ref. [24j 
in the context of EU. However, a detail study on the issue will be 
taken up elsewhere. 
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